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Recent experimental improvements on if-decay data allow for a precise extraction of the 
strangeness-changing scalar Kn form factor and the related strange scalar spectral function. On the 
basis of this scalar as well as the corresponding pseudoscalar spectral function, the strange quark 
mass is determined to be m s (2GeV) = 92±9MeV. Further taking into account chiral perturbation 
theory mass ratios, the light up and down quark masses turn out to be m u (2 GeV) = 2.7 ± 0.4 MeV 
as well as m<j(2 GeV) = 4.8 ± 0.5 MeV. As a by-product, we also find a value for the Cabibbo angle 
\V US \ ~ 0.2236(29) and the ratio of meson decay constants Fk/F^ = 1.203(16). Performing a global 
average of the strange mass by including extractions from other channels as well as lattice QCD 
results yields m s (2 GeV) = 94 ± 6 MeV. 

PACS numbers: 12.15.Ff, 14.65.Bt, 11.55.Hx 



I. INTRODUCTION 

Together with the strong coupling, quark masses are 
fundamental QCD parameters of the Standard Model 
(SM). Precise values are required as an input in many 
phenomenological applications and they can also prove 
useful in constraining unified theories which go beyond 
the SM. In the light quark sector, of special interest is the 
mass of the strange quark, as it plays an important role 
in theoretical predictions of the parameter for direct CP 
violation e'/e and SU(3) breaking in hadronic B decays. 

Present day approaches to determine quark mass val- 
ues include on a more fundamental level lattice gauge the- 
ory, whereas on a more phenomenological level, operator 
product expansion (OPE) as well as analyticity proper- 
ties can be exploited to relate the quark-gluon picture of 
QCD to hadronic observables [1], thereby also allowing 
for an extraction of quark masses. Further information 
on the light quark masses can be obtained from chiral 
perturbation theory (%PT) [2], which permits to make 
predictions for quark mass ratios. In this work, the latter 
two approaches will be considered to first determine the 
mass of the strange quark m s and from that also the up 
and down quark masses m u and mj. 

Good sensitivity to m s is achieved in the analysis of 
scalar and pseudoscalar hadronic channels, since the cor- 
responding two-point correlation functions turn out to 
be proportional to m 2 .. Both channels have already been 
studied extensively in the literature. In the scalar chan- 
nel, after the pioneering work of [3], recent analyses have 
been performed in refs. [4—11], whereas the pseudoscalar 
channel has been investigated in refs. [12, 13]. A fun- 
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damental ingredient for the scalar channel analysis is the 
scalar strange spectral function, which has been obtained 
in a series of articles [4, 14, 15] on the basis of dispersion 
relations and S-wave Kit scattering data. 

Recent experimental improvements on iC-decay data, 
to be discussed in the next section, allow for a precise 
determination of the ratio Fx / 'F„ / 'F^ '(0), and conse- 
quently for a substantial improvement of the scalar Kit 
form factor and the related strange spectral function 
which will be performed in section 3. In section 4, the 
scalar strange spectral function is employed to determine 
the strange mass m s . For comparison, we also consider 
the m s determination from the pseudoscalar channel in 
section 5. Finally, in the conclusions, our final average for 
m s will be presented and compared to recent extractions 
from lattice QCD, and from our result for m s and xPT 
mass ratios, also m u and uid are deduced. 



II. THE RATIO F k /F v /f£"(0) 

Todays knowledge of the pseudoscalar meson decay 
constants dominantly originates from the leptonic decays 
P — > Iv [16]. The ratio of the decay rates T\K — > ^(7)] 
and r[-7r — » ^(7)] is directly proportional to the square of 
Fk/Fk- Including electromagnetic radiative corrections 
according to Marciano and Sirlin [17], it reads 

rj^Mr)] = \v us \ 2 f k m k {i-xIY [i+|fM 

r[7r -> M7)] \v ud \ 2 F 2 M^ (1 - xlf [1 + 



1 



-A(7 2 ^ln 



■AC, 3, 



}<!> 



where xp = mi/Mp and ACi = Cm — C^. The con- 
stants Cip depend on the hadronic structure, and an ex- 
plicit expression for the function F{xp) can be found in 
ref. [17]. 

While extracting Fk/F^ from the ratio (1), the domi- 
nant uncertainty is due to the CKM matrix element V us . 
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Thus, it is preferable to present a value for the ratio 
\V US \F K 



\V ud \ Fir 



0.27618 (39) (27) (6) = 0.27618 (48) . (2) 



To arrive at this result, the experimental K^2 and 71^2 
decay rates have been used. For the decay rate, 

we have employed the very recent result by the KLOE 
collaboration B[K+ -> M+z^d)] = 0.6366(18) [18], and 
all other experimental inputs have been taken from the 
Review of Particle Physics [16]. AC\ has recently been 
calculated in the framework of xPT [19] with the re- 
sult ACi — Z/2 In (Mk /M n ), where the chiral coupling 
Z = 0.8±0.2 arises from the electromagnetic mass differ- 
ence of the pion. For the remaining constants, the gener- 
ous estimates AC2 = Oil and AC3 = 0±3 [17] have been 
employed. In eq. (2), the separated errors correspond to 
the uncertainty resulting from the K^ branching ratio, 
the remaining experimental inputs, being dominated by 
the K lifetime, and the radiative corrections, respectively. 
The total radiative correction in eq. (1) then amounts to 
— 3.05 (16) ^, in agreement to the result of ref. [20], ob- 
tained in a different framework. 

To extract a value for Fx / 'F„ j ' F^* '(0), we now have to 
assume inputs for \V US \ F+*(0) and \V u d\- For both quan- 
tities, we employ the results \V us \F^(0) = 0.2173(8) 
from K e z decays as well as \V u d\ = 0.9738 (3), presented 
in the very recent review [21]. (For comparison, see how- 
ever also ref. [22].) One then obtains: 



F, 



K 



1.2377 (22) (46) (5) = 1.2377 (51) , (3) 



where the different errors correspond to the uncertain- 
ties of eq. (2), \ V US \ F^ v (0), and \ V u d\, respectively. From 
eq. (3), we observe that the uncertainty on the consid- 
ered ratio is dominated by the experimental result for 
\V US \ F+*(0). The influence of \V u d\ is rather small, and 
even increasing its error by a factor of two practically 
would have no effect. 

In all previous expressions the hadronic matrix ele- 
ments F K , Fn and Ff "(0) are defined in the framework 
of pure QCD, i.e. in the limit oqed = 0. To the quoted 
level of precision, electromagnetic corrections have a cru- 
cial effect in the measured physical observables and are 
taken into account through explicit correction factors, as 
shown in eq. (1) for Fk / F^ [17, 23]. The quoted value for 
\V US \ F? n (0), has been derived from K e3 data [21], using 
the electromagnetic corrections computed in ref. [24]. 



III. THE SCALAR K it FORM FACTOR 

As the next step towards the determination of light 
quark masses, the scalar Ktt form factor F c f r ' r (i) = F (t) 
[50] will be calculated along the lines of the dispersion 
theoretic approach of refs. [14, 15], employing additional 
constraints both at long and short distances. For solving 



the system of dispersion-integral equations, in this ap- 
proach two external input parameters are required. Like 
in ref. [15], these inputs will be the value of the form fac- 
tor at t = 0, i<b(0) = F + (0), as well as its value at the 
Callan-Treiman point A KlT = M K - M%, F (A Kw ). 

For i*o(0), we employ an average over recent deter- 
minations from the lattice and effective field theory ap- 
proaches [25-30] 



F (0) = 0.972(12) . 



(4) 



which is also compatible with the original estimate by 
Leutwyler and Roos [31]. Together with the results of 
the previous section, this choice corresponds to 



IK. I = 0.2236(29) and 



K 



Fir 



= 1.203(16). (5) 



The value for \V US \ is compatible with unitarity at the 
1.2 a level, while Fk /F v is about 1 a lower than the result 
of ref. [31]. The ratio of decay constants is also in nice 
agreement with the recent lattice results [32]. 

Because the value at the origin only concerns the 
global normalisation of Fo(t), from the dispersive ap- 
proach a more precise determination of the form factor 
shape can be obtained by first concentrating on the ratio 
Fo(t)/ Fq(0). On the other hand, to a very good approx- 
imation, the second required input F (Aktt) is given by 
Fk/Ftt, which then turns into the relation 



K 



1 — = 1-2346(53) , (6) 



Fo(A K ir) = 

F Q (0) F w F o (0) ' F (0) 

where Act = — 3 • 10~ 3 has been obtained within xPT 
at order p 4 [33], and half of this value has been added 
quadratically to the uncertainty of eq. (3). 

Using this input parameter and performing an average 
over the different fits to the S-wave Ktt scattering data, as 
discussed in detail in [15], for the slope and the curvature 
of the form factor at the origin, we obtain: 



am 

F o (0) 



= 0.773(21) GeV" 



*ff(0) 
F o (0) 



1.599(52) GeV" 



(7) 

Our result for the slope can also be expressed in terms of 
the scalar Ktt squared radius, or the parameter Aq: 



6 §^ = 0.1806(49) fm 2 . 



A = M. 



F (0) 
2 ^o(0) 



F o (0) 



0.0147(4) . 



(8) 
(9) 



The value for A is compatible with our previous result 
presented in [25]. It is also in good agreement with the 
recent experimental result by KTeV [34] , but is about 3a 
lower then the corresponding ISTRA result [35], where 
however Fg'(0) was found to be compatible with zero. In 
figure 1, we display our result for the the scalar Ktt form 
factor Fo(t) as a function of y/t, while varying its value 
at the Callan-Treiman-point according to eq. (6). 
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Vt [GeV] 

FIG. 1: Result for the scalar Kit form factor Fo(t) with the 
range corresponding to a variation of the form factor at the 
Callan-Treiman-point according to eq. (6). 



IV. Ms FROM THE SCALAR CHANNEL 

The numerical analysis of the scalar channel proceeds 
in complete analogy to the previous analyses [4, 6, 11] 
where the detailed theoretical expressions can be found. 
The central equation for the extraction of the light quark 
masses from the scalar and pseudoscalar spectral func- 
tions takes the form 

so oo 

u*tfc(«)= I Pph {s)e- S / U ds+ fp th (s)e- s/u d S , (10) 



in which ^ t h(u) denotes the Borel-transformed theoreti- 
cal 2-point correlator, providing an exponential suppres- 
sion of higher-energy contributions where no experimen- 
tal information is available, and u is the so-called Borel 
variable. Expressions for ^th(u) up to 0{a 3 ) can be 
found in [4, 11] and we have also included the very recent 
result on the perturbative 0{a 4 s ) contribution [12, 36]. 

The OPE which is employed for calculating ^ t h(u) is 
valid only at sufficiently large u 3> Aq CD . As is well 
known, for the scalar and pseudoscalar channels a break- 
down of the OPE is expected to occur at a relatively large 
u around 1 GeV 2 , due to the presence of non-perturbative 
vacuum effects which go beyond the local condensate ex- 
pansion [37, 38] . Models of the correlation function based 
on instanton configurations, such as the instanton liquid 
model (ILM) [39-41], allow to penetrate into the region 
of smaller u. However, as realised e.g. in [13], at suffi- 
ciently large u « 2 GeV 2 , the ILM correction turns out 
to be rather small, and hence we shall avoid it by choos- 
ing 2 GeV 2 as a lower limit for the Borel variable, which 
furthermore also reduces the uncertainty from higher or- 
der a s corrections. 

The relation between the phenomenological spectral 
function and the strangeness-changing scalar form factor 
is given by 



where <jkp{s) is the two-particle phase space factor and 
like in ref. [4], we have also included the Krj contribu- 
tion F^'{s). A possible source of systematic uncertainty 
may be the neglect of more than 2-particle final states, 
on which we comment further below. Above the energy 
s , the spectral function is again approximated by the 
theoretical expression pth{s). 

Performing the m s analysis on the basis of eqs. (10) 
and (11), for the running mass in the MS scheme we find 



(12) 



to, (2 GeV) = 87.6 ±|;| MeV, 

at a central value s = 4.4 GeV 2 , where m s is most stable 
with respect to variations of u in the range 2 GeV 2 < 
u < 4 GeV 2 . This and all other input parameters whose 
variation induces a shift in m s of more than 1 MeV have 
been collected in table I. 



Parameter 



Value Am, [MeV] 



F (A Kn )/F {0) 


1.2346(53) 


+ 7.0 

- 5.3 


fo(0) 


0.972(12) 


+ 1.0 
- 1.1 


a s (M z ) 


0.119(2) 


-3.1 
+3.8 


0(ai) 


no O(af) 
2xO(q|) 


+ 1.8 

-1.1 


so 


3.9 - 5.5 GeV 2 


+3.2 
-2.6 



TABLE I: Values of the main input parameters and corre- 
sponding uncertainties for m s (2 GeV) in the scalar channel. 

The dominant phenomenological uncertainty on m s is 
due to the shape of the form factor Fo (t) while the value 
F o (0) only plays a minor role. On the theoretical side the 
main uncertainty arises from the variation of a s and to 
a smaller extent from unknown higher order corrections 
which are estimated by either removing or doubling the 
0(a A s ) correction. Finally, we have varied the parameter 
so in a rather generous range. Higher order corrections in 
the OPE have been included according to [4, 11]. How- 
ever, they have only small influence and a variation of 
the corresponding parameters induces errors on m s of 
less than 1 MeV in all possible cases. 



V. Ms FROM THE PSEUDOSCALAR 
CHANNEL 

In complete analogy to the scalar channel the strange 
mass m s can also be extracted from the pseudoscalar 
channel. The phenomenological spectral function has 
been modelled along the lines of ref. [13], while our anal- 
ysis parallels the one of the recent work [12]. 

From the pseudoscalar channel, at an so = 4.2 GeV 2 
for which m s in the region 2 GeV 2 < u < 4 GeV 2 is most 
stable against a variation of the Borel variable u, the 
strange quark mass is found to be 



m s (2GeV) 



97.2 ±\™ MeV. 



(13) 



(11) 



Again, in table II the input parameters and their vari- 
ations which produce a shift of m s larger than 1 MeV 
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are compiled. The most important parameters are the 
decay constants of the first two excited K resonances, 
the fC(1460) and A"(1830), which have been estimated in 
rcf. [13]. Since especially the decay constant of the sec- 
ond if (1830) resonance is not very well determined, to 
be conservative a generous range has been assumed. We 
should note, however, that employing the central input 
parameters of ref. [12], their results are exactly repro- 
duced. Furthermore, the dependence on the QCD cou- 
pling, higher order corrections and the parameter sg has 
been analysed as in the previous section. 



Parameter Value Am s [MeV] 



fk(1460) 


22.0 ± 2.6 MeV 


+ 5.2 
- 4.4 


■fjf(1830) 


lu -10 


+ 6.9 
- 4.6 


a s (Mz) 


0.119(2) 

no 0(aj) 
2xO(af) 


-4.5 
+6.6 


0(a 4 s ) 


+ 1.4 
-0.8 


so 


3.7 - 5.5 GeV 2 


+2.8 
-1.5 



TABLE II: Values of the main input parameters and cor- 
responding uncertainties for m s (2GeV) in the pseudoscalar 
channel. 



VI. CONCLUSIONS 

In the previous two sections, the strange quark mass 
has been determined on the basis of improved results 
for the strangeness-changing scalar spectral function pre- 
sented in section 3, and the resonance model [13] in the 
case of the pseudoscalar spectral function. As can be ob- 
served from eqs. (12) and (13) within the uncertainties, 
both results are in reasonable agreement, and thus, we 
are in a position to average them. Since it is difficult to 
assign a precise meaning to the theoretical uncertainties, 
we have decided to take the arithmetic mean, and as- 
signed the larger error of (12) as the total uncertainty, 
which yields our final result 

m s (2 GeV) = 92.4 ± 8.8 MeV = 92 ± 9 MeV , (14) 

providing the first determination of m s at the 10% level 
from non-lattice approaches. The fact that m s from the 
scalar channel turns out smaller than for the pseudoscalar 
channel might be attributed to the fact that only 2- 
particle intermediate states have been included. Never- 
theless, on the basis of large-Ac- arguments, the contri- 
bution from higher-multiplicity states is expected to be 
suppressed, and also the uncertainties in the pseudoscalar 
channel are so large that at present no significance can be 
attributed to this difference. Still, we plan to investigate 
this question further in the future. 

Making use of our final result on m s of eq. (14), and 
two particular quark mass ratios obtained from xPT [42], 
R = m s /m = 24.4±1.5 as well as Q 2 = (m 2 s -m 2 ) / (m 2 d - 
m u) = (22.7 ± 0.8) 2 , where rh = (m u + rrid)/2, we are 



also in a position to calculate the light up and down quark 
masses with the result: 

m u (2GeV) = 2.7 ±0.4 MeV, (15) 
m d (2GeV) = 4.8 ±0.5 MeV. (16) 

To conclude, let us compare our determinations of m s 
presented in eqs. (12) and (13) with other recent extrac- 
tions of this quantity from sum rules, and lattice QCD. 
To this end, in figure 2 the m s values obtained in this 
work, as well as the recent determinations from e + e~- 
scattering [43] and hadronic t decays [44], are displayed 
as the full circles. The most recent results from lattice 
simulations at N f = 2 + 1 [45, 46] and N f = 2 [47-49] 
quark flavours are shown as the full squares. 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

I ^ | Scalar SR [this work] 

I | Pseudoscalar SR [this work] 

I m | eV SR [Narison 05] 

I | x SR [Gamiz 05] 

I , | Lattice [MILC/HPQCD 06] 

I u | Lattice [CP-PACS/JLQCD 06] 

I , | Lattice [QCDSF/UKQCD 05] 

I u | Lattice [SPQ CD R 05] 

I B | Lattice [Alpha 05] 

I + | Average 

.. ..I .... I .... I .... I .... I .... I .... I .... I .... I .... I .... I .... I .... I .. .. 
60 70 80 90 100 110 120 130 140 
m s (2 GeV) MeV 

FIG. 2: Compilation of recent determinations of m 3 from sum 
rules (circles) and lattice QCD (squares). 

All results for m s are in good agreement - perhaps with 
the exception of the QCDSF/UKQCD value [47] which 
lies a bit high - and thus we are in a position to present a 
global average. To do this, we have calculated a weighted 
average of all numbers, taking the larger uncertainty in 
the case of unsymmetric errors, which leads to 

m s (2GeV) = 94 ±6 MeV, (17) 

where as the total uncertainty we have chosen the error 
of the single best determination of m s from the lattice. 
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